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ABSTRACT 

Context. Understanding the evolution of the dark matter halos of galaxies after they become part of a cluster is essential for under- 
standing the evolution of these satellite galaxies. 

Aims. We investigate the potential of galaxy-galaxy lensing to map the halo density profiles of galaxies in clusters. 
Methods. We propose a method that separates the weak-lensing signal of the dark-matter halos of galaxies in clusters from the weak- 
lensing signal of the cluster's main halo. Using toy cluster models as well as ray-tracing through N-body simulations of structure 
formation along with semi-analytic galaxy formation models, we test the method and assess its performance. 

Results. We show that with the proposed method, one can recover the density profiles of the cluster galaxy halos in the range 30 - 300 
kpc. Using the method, we find that weak-lensing signal of cluster member galaxies in the Millennium Simulation is well described 
by an Navarro-Frenk- White (NFW) profile. In contrast, non-singular isothermal mass distribution (like PIEMD) model provide a 
poor fit. Furthermore, we do not find evidence for a sharp truncation of the galaxy halos in the range probed by our method. Instead, 
there is an observed overall decrease of the halo mass profile of cluster member galaxies with increasing time spent in the cluster. 
This trend, as well as the presence or absence of a truncation radius, should be detectable in future weak-lensing surveys like the 
Dark Energy Survey (DES) or the Large Synoptic Survey Telescope (LSST) survey. Such surveys should also allow one to infer the 
mass-luminosity relation of cluster galaxies with our method over two decades in mass. 

Conclusions. It is possible to recover in a non-parametric way the mass profile of satellite galaxies and their dark matter halos in 
future surveys, using our proposed weak lensing method. 

Key words, gravitational lensing: weak. Galaxies: clusters: general. Galaxies: evolution, dark matter 



1. Introduction 

Astronomical observations indicate that the majority of matter in 
the Universe is of a yet unknown form that neither emits nor ab- 
sorbs Ught (e.g. [Komatsu et al.|2010{ Percival et al.|2010 Riess 



|et al . 2009). A unique way to study this so-called dark matter 
and its relation to luminous matter is provided by gravitational 
lensing ( Schneider et al.|[2006[ l. Gravitational lensing has been 
used to probe the matte r associated with galaxies (e.g. SLACS 



surve> 



Hoekstra et al. 2004, Mandelbaum et al.]l2006 Parkerj 



'et al. 2007 Mandelbaum et al. 2008 Simon et al. 2008, Tian 
et al. 2009 ), clusters (e.g. .Clowe et al. 2006, Bradac et al. 2006 



Halkola et al.|2006l|Halkola et al.|2008||Jee et al.|2009||Schirmer 



et al.|20l0|) an d the large-scale structure (e.g.|Schrabback et al. 
2010[|Fuet al. 2008 ; Benjani lneral.|2007| l 



According to the hierarchical structure formation paradigm, 
local overdensities created in the early Universe collapse into 
smaller dark matter halos, in which galaxies form. Larger halos, 
coiTesponding to galaxy clusters, form later through accretion 
and mergers. As a result, a typical galaxy cluster has a massive 
main halo of dark matter with a bright central galaxy (BCG) at 
its center Within this main halo, there are many smaller subhalos 
hosting a satellite galaxy. These were once isolated objects that 
merged with the cluster. 

An important open question is how subhalos evolve after 
they become part of a cluster This is also essential to under- 
stand the evolution of the satellite galaxies embedded in them. 
Simple analytic models assume that subhalos are just stripped 



http://www.slacs.org 



of their mass outside some tidal radius by gravitational tidal 
forces. On the other hand numerical simulations indicate that 
tidal forces also heat the subhalos causing them to expand and 
decrease their central density (e.g. Ghigna et al.||1998 Ghigna 



|et al.,2000, Hayashi et al. 2003) . In this picture, both tidal strip- 
ping and heating continually change the radial mass profiles of 
subhalos, eventually destroying them. 

Most observational studies of subhalos in clusters with grav- 
itational lensing employ parametric models for the mass distri- 
bution of the subhalos. Using parametiic models has the advan- 
tage that one can obtain useful constraints on the model parame- 
ters even with a modest amount and quality of lensing data. The 
method analyzes individual clusters and obtains the mass model 
parameters as a function of the luminosity of the satellite galaxy 



(e.g. 'Limousin et al."2007; Natar ajan et al.|2007 Halkola et alj 
2007; Suyu & Halkola 2010 ). However, a major disadvantage 
of this approach is that it relies on strong assumptions about the 
mass profiles of subhalos. 

Direct measurements of the mean tangential shear profile 
around suitably chosen samples of massive objects provide a 
more direct and less biased view of their mean mass profiles. 
This approach has been successfully applied to study the halos 
of field gal axies (e.g. [M andelbau m et al.|2008 
halos (e.g. [Sheldoii et al..,2009, pkabe et al. 



and cluster main 



2010 



Israel et al. 



|2010l. However, this approach has a drawback when used for 
cluster member galaxies and their embedding subhalos: The re- 
sulting shear signal probes the subhalo profiles only very close 
to their center, while the signal becomes dominated at larger dis- 
tances by the surrounding much more massive cluster main halo. 
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In this work, we propose a modification of the standard weak 
galaxy-galaxy lensing approach that addresses the latter prob- 
lem. The key idea is to 'calibrate out' the cluster main halo sig- 
nal by subtracting from the signal measured around the satellite 
galaxies the tangential shear signal around a specific set of cali- 
bration points. As shown in this paper, the additional calibration 
allows one to reliably probe the subhalo profiles to much larger 
radii than with the standard galaxy-galaxy lensing approach. 

We test the proposed method and assess its performance 
compared to the standard approach using simulated lensing 
fields generated from the Millennium Run (Springel et al. 2005]l. 
Furthermore, we analyze the resulting subhalos profile as pre- 
dicted by our simulations and provide forecasts for the signal-to- 
noise level expected for large upcoming and future surveys like 
the Dark Energy Survejj^DES), or the Large Synoptic Survey 
Telescop^ (LSST) survey. Moreover, we discuss how weak 
lensing can be used to study the mass loss of subhalos during 
their evolution in galaxy clusters. 

Our paper is organized as follows: In Sect. |2] we briefly 
discuss the theoretical background. In Sect. |3] we describe our 
proposed method in detail. A short overview of our simulations 
is given in Sect. |4] In Sect. |5] we test the performance of our 
method, we apply it to characterize the dark matter profiles of 
subhalos in the Millennium simulation, and we forecast signals 
and noise for upcoming surveys. In Sect. |6]we discuss the most 
relevant systematic effects for this work. The paper concludes 
with a summary in Sect. |7] 

2. Theory 

In this section, we briefly discuss the theory of gravitational lens- 
ing needed for our work. We refer the reader to the standard liter- 
ature (e.g. |Bi^E^n&jSchneider 2001 Schneider et al. 2006) 1 
for a detailed discussion. 



2.1. Weak lensing basics 

A massive foreground structure, the lens, deflects the light emit- 
ted by a distant galaxy in its background, the source, toward an 
observer. Sources at angular position = (/S^ySi) are seen by 
the observer at a possibly different image position 6 - {9\,92)- 
Differential deflection causes distortions in the images of the 
background galaxies. Locally, the distortion is quantified by the 
Jacobian 

j^ij = dPiidOj (1) 

of the lens mapping 6 t-^ /3. This Jacobian defines the conver- 
gence, 

K^l- Tr(y[)/2, (2) 



the complex shear y - yi + iyi where: 

n =(^22-^ii)/2, 

72 =(^22+^ll)/2, 



(3a) 
(3b) 



and the reduced shear g - y/(l - k). 

Defining a complex ellipticity e - e\ + iei fBartelmann' 
|& Schneider||2001) ) for the background sources, the differential 
light deflection, and thus the mass associated with the lens, can 
be inferred (statistically) from the measured shape. The relation 
between the intrinsic ellipticity and the observed one e are re- 
lated by: 



^ http : //www . darkenergysurvey . org/ 
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1 + g< 

e* +2*' 



if 1^1 < 1, 



if \g\ > 1. 



(4) 



In the following, we assume the validity of the weak lensing 
approximation i.e.: 



|/c| «: 1, g < 1, g ~ y, and e ~ e\+ y. 



(5) 



Note, however, that we test this assumption in Sect. |6.1| Under 
this approximation and assuming that the of the different 
sources are uncorrected: 



(6) 



2.2. Gaiaxy-gaiaxy iensing 



In galaxy-galaxy lensing, one measures the tangential ellipticity 
of a background galaxy image with respect to the position of a 
foreground lens galaxy. 



•cos(2(;i)ei - sin(20)e2, 



(7) 



where (p is the polar angle of the Une connecting lens and back- 
ground image. From Eq. (|6]l follows that (et) = Tt in the weak- 
lensing regime, where yt denotes the tangential component of 
the shear with respect to the lens position. 

Under the geometrically-thin lens approximation, we can de- 
fine the excess surface mass density AS(^) at projected physical 
radius ^ by: 

AS(^) = £(^)-E(a (8) 

where £(^) denotes the average surface mass density within a 
circle of radius ^, and S(f) is the average surface mass density 
on that circle. 

The tangential shear yt(^) averaged on a circle with physical 
projected radius ^ around a lens galaxy is related to the excess 
surface mass density AX(^) by ( Schneider|2006 i: 



AE(^) = rt(^)2crit. 



(9) 



where the critical surface mass density Sdit is a function of the 
angular diameter distance to the lens D^, to the source D^, and 
the distance between lens and source Da,: 



^rrir — 



(10) 



An important aspect of galaxy-galaxy lensing is to com- 
bine the different measurements from different foreground- 
background galaxy pairs to increase the signal-to-noise. Each 
lens-background galaxy pair / with projected separation ^ at the 
lens, tangential ellipticity et,; of the background galaxy image, 
and critical surface mass density Sent./ provides an estimate 



A2;(^) = Scrit.jet,,- 



(11) 



for the excess surface mass density AE(^) of the lens. For a sam- 
ple of lens and background galaxies, one can compute the final 
estimate AZ(^) by a weighted mean: 



(12) 



where the sum runs over all pairs formed by a background image 
and a lens galaxy. If one assumes that each lens galaxy has the 
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same mass profile and each shear estimate from a background 
galaxy image carries the same uncertainty, the optimal weights 
are given by 

^ciL-fc,Zl) 

Wi = — '-^ . (13) 

The weights make our estimator sensitive to the redshift distri- 
bution of lenses and sources. These weights are a simplified ver- 
sion of those used, e.g, by Mandelbaum et al. (20081, since we 



neglect redshift errors and redshift-dependent errors in the ellip- 
ticity estimation. 

In this work we also consider simulated pixelized surface 
mass density maps. In those cases, we estimate directly: 

k 

around each mass concentration. The sum runs over all the lenses 
in the sample. For these weights we use the number of pixels in 
the annulus defined to compute 



(15) 



2.3. Analytic models for halo profiles 

In the analysis of the simulated excess surface mass profiles, we 
consider several analytic models for the mass distribution inside 
galaxy halos. 



2.3.1. Navarro-Frenk-White profiles 

On average, the three-dimensional density profiles of the virial- 
ized inner regions of dark-matter halos in simulations are well fit 
by Navarro-Frenk-White (NFW) profiles ( [Navarro et al.|1997| l: 

5c Pm 



p(r) = 



(r/r,)(l+r/r,)2' 



(16) 



where 5c is a characteristic density, is the mean matter density 
of the universe at the halo redshift, and is the halo scale radius. 

We characterize the spatial extent of NFW halos by r2oo, i-C- 
the radius in which the average density is 200 times the mean 
cosmic matter density p^. A commonly used parameter for NFW 
profiles is the concentration c - r20{)/?"s- 

The mass M200 within r2oo is then given by 



M20() = AndcPmr^ 



ln(l +c)- 



1 



1 +c 



(17) 



Using the parameter x - ^/r^ and the characteristic surface 
mass density - rs6cPm, the projected surface mass density 
Snfw of the NFW profile at a projected radius ^ reads E(^) = Es ■ 
fi^lfrms) and the mean projected mass density £(^) - "Es-gi^/r^) 
where f(x) and g(x) are defined in Wright & Brainerd ( 2000[ l, 
(see also |Bartelmann| 1 996| ). 

2.3.2. Truncated Navarro-Frenk-White profiles 

If one integrates the mass of an NFW profile up to an infinite ra- 
dius, the total mass diverges unless one introduces an additional 
truncation. The truncated NFW profile we consider here was de- 
rived by Baltz et al. ( 2009| l. Its density profile reads: 



ScPn 



(r/r,)(l+r/r,)2[l+(rK.)2]' 



(18) 



where denotes the truncation radius. Defining x - ^/r^, t 
r^jr^, and the function 



F{x) 



ln(l/x+ yjllx^ - 1] 



Vl -x2 



1 



arccos(l/x) 
, Vx2- 1 



for X < 1 , 
for X = 1 , 

for jc > 1 , 



(19) 



one has 



2(x) 



(t2 -I- 1)2 

n 



1 



Vt2 + x2 T Vt2 



(1 -F(x)) + 2F{x) 

log! 



r2-l 



+ x^ 



V72 



+ X^ + Tl 



(20) 



The mean density within a radius x is: 



E(x) = 



42s 



x2 (t2-h1)2 

+ \+ 2{x^ - 1) 



F{x)+T7I+{t^ - l)logT 



+ x-^ 



1 



-n + 



log 



+ X^ + Tl 



(21) 



2.3.3. Truncated Non-singular Isothermal profiles 

This profile (in the following TNSI) is a truncated and spheri- 
cal version of the Pseudo-Isothermal Elliptical Mass Distribution 
(PIEMD) profile by |Kassiola & Kovner] ( |1993| l. Its three dimen- 
sional distribution reads: 



p{r) 



pa 



(l+'-Vdre)(l+'-Vdt) 



(22) 



with Tcore denoting the profile core radius, and rcut denotes a cut- 
off radius. The cutoff radius r^^i plays a role similar to the trun- 
cation radius rn of the truncated NFW. The surface mass density 
distribution reads: 



2:(^) = 



1 



1 



(23) 



By integrating, one obtains the mean surface mass density 



2po r; 



2 2 
r 7T 



^ cut ^ core) 



/'■cut + 



('"ci 



(24) 



3. Galaxy-galaxy lensing for subhalos 

Fig.[T]shows a schematic view of a galaxy cluster consisting of a 
cluster main halo and a subhalo. When galaxy-galaxy lensing in 
its standard form is applied to probe the mass distribution around 
satellite galaxies (subhalos), the resulting lensing signal has a 
strong contribution from the cluster main halo. 

One can use an analytic model of the main halo to take 
its contribution into account. However, this relies on strong as- 
sumptions about the main halo mass profile. We propose a non- 
parametric scheme to estimate the main halo contamination and 
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Fig. 1: Sketch of a cluster with the main halo and a subhalo with a 
cluster galaxy at its center Also shown is the calibration point used 
to estimate the main halo contribution to the lensing signal around the 
cluster galaxy and the same distance dus from the cluster center. 
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subtract it from the signal. As illustrated in Fig. |2] a calibration 
point can be used to estimate the contamination from the main 
halo. The calibration point is located at the same distance d^s 
from the main halo center as the subhalo center, but on the op- 
posite side. If the cluster is point-symmetric around its center, 
we can subtract the signal around the calibration point from the 
signal around the subhalo, and remove the contamination from 
the main halo. In the following, we will call this method galaxy- 
galaxy lensing with calibration (GGLC). 

Note, however, that the calibration does not work for very 
large projected radii ^. For ^ > 2xiiM-s the method fails since the 
calibration signal includes the subhalo signal itself. Moreover, 
as shown in Fig. [2j the signals around the subhalo and the cali- 
bration point are very steep for ^ around dM-s ■ Thus the recon- 
structed subhalo signal there will be very susceptible to errors 
introduced by shot noise or the intrinsic ellipticity of the back- 
ground galaxy. This makes it impossible to reconstruct the sub- 
halo profile for ^ > c/m-s in practice. 

Since the point-symmetry of any cluster sample is not per- 
fect, our calibration is not only limited from the theoretical point 
of view. We discuss this in detail in Sect. 15. II 

The excess surface mass density AEsubh(^) measured around 
a sample of lens galaxies can be thought to consist of (at least) 
two parts. The first contribution is due the (sub)halos of the 
lenses we are targeting. This part usually dominates the signal 
close to the lens galaxy. There, A2subh(^) can be interpreted as 
the excess surface mass density profile of an average lens galaxy 
halo. The second contribution is due to the (sub)halos of other 
galaxies, and the clusters' main halo. This contribution domi- 
nates on larger scales. There, AEsubh(^) cannot be interpreted di- 
rectly as the average lens galaxy halo profile. 

For a reaUstic study we need a large sample of clusters, and 
we do not need each cluster to be perfectly point symmetric. 
We only need the cluster on average to be sufficiently point- 
symmetric. In the following the lensing results we present are 
given by 



10"' lO"" 10' 

Physical radius [Mpc] 

Fig. 2: Galaxy-galaxy lensing signals of cluster main halo and subhalo 
(assuming NFW profiles). The subhalo is at 0.6 Mpc from the cluster 
center. Top panel: the subhalo around the subhalo center (thin solid red 
line), the main halo around the main halo center (thick solid black line), 
and the main halo around the calibration point (dashed/dotted line for 
positive/negative part). Bottom panel: the subhalo around the subhalo 
center (thin dotted green line), the subhalo plus main halo around the 
subhalo center (black dashed line), the subhalo plus main halo around 
the main halo center (cyan dotted line), the subhalo plus main halo 
around the calibration point (red solid line), and the reconstructed sub- 
halo (blue thick solid line). 



A^i^) = AEsubh - A2eal(^). 



(25) 



4. Simulations 



In order to test our proposed galaxy-galaxy lensing method and 
provide forecasts for the expected signal, we use simulated lens- 
ing fields. The fields are obtained by ray-tracing through the 
Millennium Simulation in combination with mock galaxy cat- 
alogs from semi-analytic galaxy formation models. 

4.1. The Millennium Simulation 



The Millennium Simulation (MS) (I Springel et al.||2005"] l is a 
large A^-body simulation of cosmic structure formation in a 
flat ACDM universe. It uses 2160^ dark-matter particles of 
8.6 X 10**/;"' Mq in a periodic simulation volume of 500/1"' Mpc 
side length. The assumed cosmological parameters are: a Hubble 
constant h = 0.73 (in units of 100km s"' Mpc"'), a matter den- 
sity - 0.25 and a cosmological constant energy density 
Q.A = 0.75 (in units of the critical density), and a normalization 
of the matter power spectrum erg = 0.9. 

An important aspect of the simulation is how structure is 
defined. All particles which are closer than a certain threshold 
value are linked. Independent groups of particles are grouped 
with this Friends of Friends (FOF) algorithm. These FOF groups 
are considered main halos in the simulation. Substructures 
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Fig. 3: Projected mass maps of a cluster in the MS (left panel) and the synthetic cluster created from it (right panel). 



within these main halos are identified by the SUBFIND algo- 
rithm ( jSpringel et al.|[2001| l: Each FOF group is subjected to 
a gravitational binding analysis, and every group of at least 20 
gravitationally bound particles forms a subhalo. The algorithm 
also estimates a mass Msubf for each subhalo that equals the 
total mass of the bound particlesj^ 

We define a cluster as a FOF group whose main halo 
has a mass according to the SUBFIND algorithm larger than 
10''* M0//1. We use this criterion to certify that the selected FOF 
groups have a large fraction of their mass gravitationally bound 
and to avoid non-virialized objects. 

The MS itself does not take into account baryonic physics. 
Galaxy formation within the evolving matter distribution of the 
simulation can however be followed by semi-analytic galaxy for- 
mation models (see, e.g., |Springel et al. 2001; |De Lucia et al. 
2004]!. In this wo rk we use the semi-analytic catalogues by De 
Lucia & Blaizot ( |2007[ l to infer the properties of galaxies in the 
MS. The model used by them assumes that at the center of ev- 
ery main halo and subhalo there is a galaxy whose evolution is 
determined by the merging history of the host dark-matter halo 
and gas-physical processes. These processes are approximated 
by a set of empirical equations with parameters tuned to make 
the resulting galaxy population consistent with many observed 
galaxy properties. 

Furthermore, the semi-analytic model allows galaxies whose 
subhalos has been dissolved during a halo merger event onto a 
larger halo to survive for some time before merging onto the 
central galaxy of the main halo. The properties and merging 
timescales of these orphan galaxies are not known very accu- 
rately (from observations or theory), but their contribution to the 
lensing signal can be significant. In our primary analysis we ne- 
glect them. In Sect. |6.2| we discuss how our results are changed 
if these galaxies are included in the analysis. 

4.2. Lensing simulations 

The lensing effects caused by the matter structures in the MS 
are calculated by the multiple-lens-plane ray-tracing algorithm 



described in Hilbert et al. (2009 1. The ray-tracing is used to cre- 
ate 128 simulated fields of view of 2 x 2deg^. Convergence and 
shear in these fields are stored on regular grids of 4096^ pix- 
els in the image plane for multiple source redshifts. In addition, 
the lensed image properties of all semi-analytic galaxies in each 



Note, that there is no unique way of estimating the mass of a sub- 
halo, and that even the same algorithm gives different values depending 
on the subhalo position inside the main halo ( [Muldrew et al.|2010^ . 



field with stellar masses > 10^/;"' M© are computed. The lens 
galaxy population for the galaxy-galaxy lensing simulations are 
selected from these lensed mock galaxy catalogs. 

The ray-tracing simulations are also used to compute pix- 
elized maps of the projected surface mass density "Lifi) of thin 
redshift slices of the simulated lensing light cones. We use these 
maps first to create simplified mock clusters for performance 
tests of the calibration method; and second to make a prelimi- 
nary analysis of our signals. 

The source galaxy populations for the galaxy-galaxy lensing 
simulations are created by randomly drawing image positions 
within the simulated fields (assuming a uniform spatial distribu- 
tion in the image plane). The redshifts of the source galaxies are 
drawn from a distribution ( |Brainerd et al.|1996) 



p(z) - z" exp 



{-(z/zof) 



(26) 



where a - 2, /3 - 1.5 and zo is a function of the median redshift 
of the survey. We then identify the lens plane that is closest in 
redshift to each galaxy and interpolate onto it the shear and mag- 
nification from the ray-tracing. We add an intrinsic ellipticity to 
the shear in order to mimic real galaxy images (Eq.|4|. The mod- 
ulus of the intrinsic ellipticity is drawn from a Gaussian distri- 



bution with zero mean and standard deviation cr^^.^. In Sect. |5. 3. 2 



we discuss in more detail our simulated lensing catalogues. 
4.3. Mock cluster mass maps 

To test our calibration scheme, we use a sample of mock clus- 
ters. This allows us to know the 'true' mass profiles of main halo 
and subhalos, and thus to quantify how well our method can re- 
trieve these profiles. Since we want to assess the performance 
of the method in situations as realistic as possible, we use the 
cluster sample we have from ray-tracing through the Millennium 
Simulation as a starting point. For each cluster in our cluster 
sample we create mock clusters at various levels of simplifica- 
tion consisting of an elliptical, or non-elliptical main halo, or no 
main halo at all. 

The main halo of each cluster in the ray-tracing catalog is 
characterized first by fitting an elliptical truncated NFW pro- 
file to the projected mass profile of the cluster]^ The resulting 
mass distribution 2eiiip(jS) is point-symmetric. To create non- 
symmetric main-halos, the original projected mass T-ifi) of the 



' Projected NFW profiles have been shown to be a good fit for pro- 
jected cluster halo profiles of the Millennium Simulation ( [Hilbert &] 
White|2010^ . 
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cluster is first divided by the best-fitting elliptical truncated NFW 
profile SeiiipOS). We then calculate the coefficients 



Sub-halos truncation radiur r^^ = 2 r. 



31„ - ^ 



p 2ellip(^p) 



COS 



Sin 



[n<f>iPp)] . 
[ncpifip)] , 



(27) 



where the index p runs over all pixels in the considered patch. 
The angle 0(j8^) is the polar angle of the pixel from the center of 
the cluster The synthetic non-elliptical mass distribution is then 



l-f^{cos [ncf>(fip)]+sm[n<p(J3p)]} 



(28) 



The mock subhalo profiles are based on the projected mass 
profiles of subhalos in the MS. We group subhalos in our lens- 
ing simulations into bins of subhalo mass Mjubf- Using the pro- 
posed calibration method on the projected mass maps, we calcu- 
late the bin-averaged excess surface density AE(^) as a function 
of radius ^ for each bin. We fit NFW profiles to the resulting 
profiles in the range 0.03 Mpc < ^ < 0.2 Mpc. The relations be- 
tween the subhalo mass Msubf and the best-fitting parameters Ss 
and of the NFW profiles are then fit by power laws, yielding 



= 1.1x10 and 

Mpc \ Mo j 



^^.1.92x10-^(^1 
MoPc-2 \ Mo / 



0.24 



(29) 
(30) 



The mass of the NFW profile without truncation diverges if we 
integrate up to an infinite radius. Therefore, we truncate our 
halos to make them more realistic. We try different truncation 
schemes, keeping the former relations for and Ej. We experi- 
mented with subhalos without truncation, significantly truncated 
(rtr = 2 Ts) and with a truncation above the maximum radius for 
which we can recover the subhalo lensing signal (rn - 6.66 r^) 
in the least massive subhalos. 

Each synthetic cluster is constructed by combining the pre- 
viously described mock main halo and the mock subhalo pro- 
files, whereby the centers of the main halo and the subhalos are 
used as centers for the mock halos. The mock clusters are then 
put onto grids having the same resolution as the projected mass 
maps of the original clusters. Thanks to our procedure, the mock 
clusters properly take into account 

- the masses and sizes of the main halo and the subhalos, 

- the spatial distribution of the subhalos within the cluster, and 

- the correlation between the main halo shape and the position 
of the halo center and the subhalos. 

An example cluster is shown in Fig. |3] 
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Fig. 4: Ratio of the measured AS(f ) to the one we input, using GGLC 
and standard GGL. The dotted black line, where the ratio is one, is a 
visual reference. 
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5. Results 

5. ^ . Performance tests 

The method we proposed in Sect.|3](GGLC) allows theoretically 
to clean the galaxy-galaxy lensing signal of a subhalo sample 
from the main halo contamination. Here, we investigate the per- 
formance of the method using the synthetic cluster sample. We 



Fig. 5: Results from mock clusters. Ratio of the measured AE(<f) to 
the one we input using GGLC, versus the projected physical radius f . 
Top figure: non-truncated NFWs. Bottom figure: = 2r^. We plot 
only the results for a selection of mass bins, the mass ranges are given. 
Black dotted-dashed lines: no main halo was added. Red dashed lines: 
clusters with an elliptical main halo. Blue solid lines: clusters with a 
non-elliptical main halo. The dotted black line, where the ratio is one, 
is a visual reference. 
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quantify the range where we can recover the signal below a cer- 
tain bias and compare it to the reliably range when using the 
'standard' GGL. 

We divide the whole mock subhalo sample into several bins 
according to Msubf as listed in Table [T] For each bin, we calcu- 
late the stacked excess surface mass density profiles AS(^) as a 
function of projected radius ^ one would obtain if the subhalos 
were completely isolated. We then compare this 'true' profile to 
the signal measured from our synthetic clusters. 

In Fig.|4j we present the difference between GGL and GGLC 
for one of the truncation schemes. The different truncation 
schemes offer the same qualitative results. The bias with respect 
to the true value is reduced in all cases. For the most massive 
subhalos, comparable in mass to the main halo, the improvement 
is reduced. On the other hand, for the least massive subhalos the 
range in radius ^ with a small enough bias doubles. This is a 
significant improvement for those subhalos which are more sub- 
jected to tidal fields and therefore better tracers of the effects that 
we want to study. 

Fig. |5] shows the ratio of the measured and true signal for 
samples of subhalos in the clusters with truncated and non- 
truncated NFW subhalo profiles and various models for the 
mass distribution of the cluster main halo. For each mass bin 
in each panel, we present the results from different situations. 
We present AE(^) for the following cases: coming only from the 
subhalos, i.e. the main halo was not added; for an elliptical main 
halo plus the subhalo population; and from a realistic cluster, 
using a non-elliptical main halo. 

The distribution of the sub-halos is not perfectly point- 
symmetric with respect the cluster center. When we measure 
around a particular subhalo, we "see" the surrounding ones and 
the contribution to our measurement does not average out. This 
also shows a correlation of the subhalo positions independent of 
the mass distribution of the main halo. 

A completely elliptical halo can always be calibrated prop- 
erly. The observed difference with respect to the signal with sub- 
halos only, is due to the resolution of our maps (which we also 
use for the synthetic cluster sample), combined with the sensitiv- 
ity of the calibration for radii larger than ^/m-s (Sect.[3]l. Internal 
checks were able to verify the source of this discrepancy. 

As a compromise between a small bias and a sufficiently 
large range in ^, we measure our profiles for those radii where 
the bias is less than a 10% from the input model. We will only be 
able to detect a truncation radius if this scale is well inside the 
range we can measure. The difference between the truncation 
models does not substantially change the region that we even- 
tually define as suitable to derive the mass profile. Since as we 
show later on, our simulations do not present a detectable trun- 
cation of the mass density profile, we focus in the following on 
subhalos without truncation. 

With GGLC, the measured subhalo profile recovers the true 
profile within 10% accuracy up to a radius of ~ 200 kpc for 
less massive subhalos (with Msubf ~ 10"'^/i"' Mq and rs ~ 
0.03 Mpc), and up to ~ 500 kpc for more massive subhalos (with 
Msubf ~ lO'^^/i"' Mq). We can also see that on scales below 
40 kpc, the measured signal is affected by the limited resolution 
of the synthetic mass maps. 

The main halo's contribution to the bias in AS(^) is more 
important for subhalos with a small separation from the cluster 
center cIm-s- In fact, including subhalos with a too small liM-s, 
may decrease the quality of the measured AZ(^). Having this in 
mind, we consider different minimum values for c/m-s • 

The results are shown in Fig. |6] Here we present again the 
ratio of the measured A2(^) to the input value versus physical 



Table 1: Binning according to Msubf of all subhalos in clusters for 
z < 0.9: Mass ranges, number of subhalos, and mean masses of each 
subhalo bin used for assessing the performance of GGLC. Note that 
mass ranges are given in /i^'Mq, while mean masses are given in M©. 



Bin 

r 

2 
3 
4 
5 
6 
7 



Mass range [Mo/h] 



[10"^' 

[1011.66 

[10'2 

[1012.33 
[1012.66 

[10'3 

[1013.33 
[1013.66 

[10'" 



10'2 ) 

1012.33) 
1012.66) 

10" ) 
1013.33) 

1013.66) 

10'"* ) 
1014.33) 



#Subhalos 

3T831 
16025 
8273 
3834 
1816 

810 

385 

162 

74 



Msubf [Mg] 
4.27 X 10" 

9.14 X 10" 
1.97 X 10'^ 
4.23 X 10'^ 

9.15 X 10'- 
1.96 X 10'^^ 
4.27 X 10'^ 
8.92 X 10'^ 
1.88 X 10"' 



Table 2: Range in radius ^ for which the bias is below 10%, for each 
subhalo bin according to Msubf in clusters for z < 0.9. The subhalos 
have a separation from the main halos center rfn-s > 0.5 Mpc. The 
mass binning can be found in Tab.[T] 



Msubf bin 


^ range 


1 


[0.043:0.18] 


2 


[0.043:0.18] 


3 


[0.043:0.18] 


4 


[0.043:0.27] 


5 


[0.043:0.27] 


6 


[0.043:041] 


7 


[0.043:041] 


8 


[0.043:0.5 ] 


9 


[0.043:0.5 ] 



radius ^. This time we compare different selections of subha- 
los. The signal improves if we discard subhalos with c/m-s < 
0.5 Mpc. The most massive subhalos are more likely to be recent 
mergers and lie far away from the cluster center, so the effect is 
small in this case (lowest panel). 

Increasing the minimal cIm-s, depletes the sample and lim- 
its our result to subhalos in the outskirts of the cluster. For 
'^M-s < 1 Mpc the range where the bias is below 10% increases, 
but only at the expense of depleting our sample. For the less 
massive halos, the depletion of the sample worsens the signal. 
For the most massive subhalos, a larger minimal c/m-s improves 
the results, and seems to be a better choice. Nevertheless, in this 
work we shall discard subhalos with d^s < 0.5 Mpc, regardless 
of subhalo type for simplicity. We present the range in ^ defined 
with our mock clusters for all Msubf mass bins in Tab.|2] 

5.2. Subhalo-galaxy lensing in the Millennium simulation 

5.2.1 . Subhalo mass profiles for different Msubf 

We now analyze the mass profiles of subhalos in the MS as ob- 
tained by applying the GGLC to the projected mass maps from 
the simulation. The subhalo mass ranges, numbers, and mean 
masses of the subhalo samples used are listed in Table [3 For 
the analysis, we only consider subhalos in clusters with redshift 
z < 0.9 and with c/m-s < 0.5 Mpc. The resulting average excess 
surface mass density profiles are shown in Fig.|7] 

In order to characterize the subhalo mass profiles we try to 
find the best fitting model. We consider an NFW profile, a trun- 
cated NFW profile and a TNSI profile. For the NFW model, we 
choose the characteristic surface mass density and the scale 
radius as free parameters. For the truncated NFW model, we take 
the truncation radius as additional parameter For the TNSI 
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Fig. 6: Results from mock clusters, ratio of the measured AZ(f ) to the 
one we input, versus the projected physical radius f . Only the results 
for a selection of mass bins, for the case with non-truncated NFWs are 
presented. The mass ranges are given. Blue dashed lines: whole subhalo 



sample. Solid red lines: subhalos with rf^ 



0.5 Mpc. Black dotted 



lines: subhalos with d^s > 1 Mpc. The dotted black line, where the 
ratio is one, is a visual reference. 

model, we fix the core radius rcore = 10""* Mpc and only fit the 
cut radius rcut and the central density po. This allows us to ana- 
lyze the cho ices made by Limousin et al. ( 2007| ) and |Natarajan 
|^aL] ( [2007l l. 

We conduct a Bayesian analysis to get the best fitting param- 
eters for each model in each mass bin. We assume a Gaussian 
likelihood in each case. The corresponding covariance matrices 
are computed from the data using a bootstrap algorithm. We 
construct each bootstrap realization by randomly drawing sub- 
halos with replacements. We do so for each mass bin, producing 
10000 realizations in each case. 

Table 3: Binning of all subhalos in clusters for z < 0.9 according 
to MsuBF- All subhalos with c/m-s > 0.5 Mpc. Note the samples are 
slightly different to those in Tab.[T] The ranges are given in Mg/h while 
the mean mass are in M^. 



Bin 


Range [M^/h] 


# Subhalos 


MsuBF [Mq] 


1 


[1011.33 


1011.66) 


23229 


4.27 X 10" 


2 


[1011.66 


10'2 ) 


11821 


9.14 X 10" 


3 


[10'2 


1012.33) 


6151 


1.98 X 10'^ 


4 


[1012.33 


1012.66) 


2837 


4.25 X 10'^ 


5 


[1012.66 


10" ) 


1476 


9.22 X 10'^ 


6 


[10'3 


1013.33) 


728 


1.97 X 10'^ 


7 


[1013.33 


1013.66) 


363 


4.30 X 10'-'' 


8 


[1013.66 


10" ) 


155 


8.79 X 10'^ 


9 


[10'" 


1014.33) 


73 


1.89 X 10'" 



We use the Bayesian evidence to rank the different models. 
The evidence was computed with a nested sampling algorithm 
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Fig. 7: The AE(f ) profile measured from the projected mass maps 'L(J3), 
for the different mass bins. Thick lines show the final range used for the 
analysis. The range in ^ was defined through the test in Sect. |5.1| In 
dashed black we show the average signal for the host cluster as a visual 
reference. 

Table 4: Comparison of the logarithm of the evidences for the three 
models fitted to the AX signals. In bold we emphasize the best model in 
each mass bin. 



Bin 


NFW 


trunc. NFW 


TNSI 


1 


-89.71±0.05 


-89.57±0.08 


-157±1 


2 


-88.31±0.09 


-88.54±0.09 


-285±1 


3 


-88.38±0.09 


-89.38±0.12 


-472±1 


4 


-85.90±0.07 


-87.72±0.13 


-710±3 


5 


-98.32+0.07 


-101.68±0.12 


-876±3 


6 


-93.38±0.05 


-97.04±0.17 


-576±1 


7 


-99.20±0.06 


-102.63±0.21 


-809±2 


8 


-88.09±0.06 


-92.80±0.18 


-528±1 


9 


-85.62±0.04 


-89.42±0.14 


-445±1 



( |Skilling|[2004) l. As shown by the results in Table |4] a TNSI is 
not a suitable model for the subhalos in the MS. Our simulation 
does not contain baryons, which are expected to change the mass 
profiles at the innermost regions. Nevertheless, our results raises 
questions about the choice of the PIEMD model for lensing ob- 
servations of cluster galaxies. The best model (except for one 
mass bin, which we consider to be due to statistical fluctuations) 
is a simple NFW profile. The values for the truncation radii ob- 
tained by the fits were in both cases (TNSI and truncated NFW) 
larger or comparable to the virial radius of the average cluster, 
and therefore consistent with a lack of a detection. 

The errors presented are the standard error from 10 realiza- 
tions of the measurement except for the leftmost model where 
we use a raw estimate coming from the algorithm. The analysis 
in this case did not require a better accuracy. We write in bold 
the maximum in each case, which belongs to the best model. 

From our simulations, it is not possible to detect a trunca- 
tion radius. A truncation, if present at all, happens at ^ larger 
than 0.2 Mpc, where our measurements become unreliable due 
to bias. This is up to four times larger than the truncation radii 
infe rred by ,Limousin et al. (20071, Natarajan et al. ( 2007| l and 
Halkola et al. '{ 2Q07\ from lensing data. 

We take the mean of the posterior from the nested sampling 
algorithm as the best-fitting values for the model parameters, and 
the posterior standard deviation as parameter errors. The best- 
fitting values for the concentration c and M200 assuming an NFW 
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Fig. 8: Parameterization of the subhalo profiles using NFW profiles. 
Bottom panel: concentration concentration as a function of Msubf- Top 
panel M200 as a function of Msubf- 

Table 5: Binning of all subhalos according to Msubf in clusters for 
z < 0.9. All subhalos with ^m-s > 0.5 Mpc. The mass binning can be 
found in Tab.|3] 



M. 



'suBF bin 


Fitted M200 [Mq] 


1 


(9.56 


±0.16)x 10" 


2 


(2.25 


±0.03)x 10'2 


3 


(4.51 


± 0.09) X 10'^ 


4 


(1.058 


± 0.016) X 10" 


5 


(2.16 


± 0.04) X 10" 


6 


(4.31 


±0.11)x 10" 


7 


(8.10 


±0.16)x 10" 


8 


(1.56 


± 0.07) X 10'" 


9 


(2.24 


±0.09)x lO'* 



profile are presented in Fig. 8] The concentration c slightly de- 
creases with increasing subhalo mass Msubf- The mass M200 
shows a tight and essentially linear relation with the subhalo 
mass Msubf- Thus, M200 is a good observable to infer the sub- 
halo mass Msubf- Note, however, that the mass M200 is always a 
few times larger than Msubf- This indicates that the subhalos as 
defined by SUB FIND do not fully extend to the radius r2oo- 

The use of M200 as a direct mass estimator on subhalos 
should be avoided. Otherwise a large fraction of the mass not 
bound to the subhalo (e.g. nearby subhalos) is attribute to it. We 



Table 6: Best fitting parameters for the relation Msubf(M2oo)/Mo = 
A (M2oo/Mpc)" coming from the data displayed in Fig.|9]in the bottom 
right panel. We present the mean and the standard error of the posterior 
of the parameters. 



Data 


n 


logio(A) 


Mass bin. 


1.059±0.013 


-1.13 ±0.17 


r band bin. 


1.20±0.04 


-2.5 ± 0.4 


i band bin. 


1.18±0.03 


-2.5 ± 0.4 



use M200 only as an observable from which we can derive the 
gravitationally bound mass Msubf. 

In observations, one does not know the radial range at which 
the GGLC signal can be reliably measured beforehand, since the 
range depends on the subhalo masses to be inferred from the 
lensing signal. Here we use the following method to resolve this 
problem. We infer M2()o from the measured GGLC signal in a 
preliminary radial range, i.e. up to 1 Mpc. Based on the results 
in Tab. [5] we use a smaller more adequate range in ^ to measure 
M200 again and iterate until obtaining consistent results. 

5.2.2. Subhalo mass estimation for magnitude-selected 
subhalo samples 

So far, we have only studied mass-selected subhalo samples. We 
now explore how well one can estimate subhalo masses Msubf 
from stacked subhalo samples that are selected according to the 
magnitudes of the hosted galaxies, which are more amenable to 
observation. 

Figure. |9]shows the mass M200 obtained from the GGLC ver- 
sus the true mean subhalo mass Msubf for sub-halos binned ac- 
cording to Msubf, observer-frame absolute r-band magnitude, 
and observer-frame absolute /-band magnitude. We define 14 
bins in luminosity, with a bin width of 0.33 magnitudes, between 
the values -19 and -23.66 for the luminosity or both bands. The 
relation between the mean inferred mass M200 and mean sub- 
halo mass Msubf is roughly linear, but noticeably depends on 
the quantity used for binning. We characterize the relation be- 
tween M200 and Msubf by a power law. 



Msubf 
Mo \M^j- 



(31) 



with best-fitting amplitudes A and exponents n listed in Table [6] 
The mass-binned and magnitude-binned subhalo samples 
differ, in particular in the distribution of subhalo masses and 
mass profiles. As a consequence, even for the same mean sub- 
halo mass Msubf, the mean profiles may differ substantially. 
While the averaging to obtain the mean Msubf and the mean 
GGLC for each bin is a linear process, the computation of the 
best fitting parameters and M200 is not. Thus, we may not be 
able to estimate Msubf with high precision through M200, but 
we can derive useful constraints on it. 

5.2.3. Subhalo evolution 

In order to observe how subhalos evolve inside the clusters we 
considered two new parameters, the infall mass Mjnt and infall 
time (both computed by [Wang et al.|20(j6l l. The first is the mass 
of the main halo that evolved into the current subhalo at the last 
epoch when it was a central dominant object, i.e. the mass just 
before it was accreted by a larger structure. The infall time is the 
moment at which this accretion was detected and allows us to 
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Fig. 9: Fitted M200 as a function of average Msubf for different bin- 
nings using the projected mass maps S()S). Black solid line: binning 
according to Msubf- Dotted red line: binning according to r band ab- 
solute observer frame luminosity. Dashed green line: binning according 
to ( band absolute observer frame luminosity. 

Table 7: Properties of subhalo samples used to study the time evolu- 
tion of subhalos: symbol and line type used in plots, subhalo numbers, 
average ^m-s in Mpc, mean subhalo mass Msubf^ and mean inferred 
M200- Only subhalos in clusters with redshift z < 0.9 and with distance 
to the main halo center d < 0.5 Mpc. 



Bin 


# Subii. 


<^M-S 


AfsUBF [Mq] 


M200 [Mq] 


X — 


31408 


1.149 


2.77 X 10" 


(6.41 ± 


0.13) X 10" 


+ — 


14927 


1.132 


1.68 X 10" 


(4.95 ± 


0.17) X 10" 


• — 


10482 


1.175 


(1.236 X 10")* 


(3.27 ± 


0.17) X 10" 


X- - 


11487 


1.169 


7.38 X 10" 


(1.85 ± 


0.03) X 10'- 


-I-- - 


6391 


1.136 


4.29 X 10" 


(1.23 ± 


0.04) X 10'^ 


• - 


4410 


1.143 


2.9 X 10" 


(7.5 ± 


0.3) X 10" 


X ■ ■ ■ 


10792 


1.260 


8.1 X 10'^ 


(1.28 ± 


0.03) X 10'' 


H 


4782 


1.151 


2.1 X 10'^ 


(6.09 ± 


0.17) X 10'- 


«... 


2331 


1.150 


1.04 X 10'^ 


(2.98 ± 


0.11) X 10'^ 



'The value after eliminating two misidentified subhalos. 

compute the subhalo age, i.e. the time spent inside the cluster at 
the time of its observation. 

We divide our sample into three classes according to Minf, 
and each of these subhalo classes into three sub-classes accord- 
ing to their age. In order to define the classification, we maxi- 
mized the signal-to-noise defined as the average ratio (over ^), 
between the measured AS(^) and its sample variance (we also 
avoided too small subhalo samples). The resulting ranges for 
classes and sub-classes are listed in Tab.|7] The resulting GGLC 
signals are shown in Fig. [TO] 

As shown in Table|7] the mean subhalo masses decrease with 
increasing subhalo age independent of the considered infall mass 
range. This is clear evidence for subhalos losing mass inside 
clusters|^The mass loss is also seen in the subhalo GGLC signals 



shown in Fig. 10 which also decrease in amplitude at smaller 
radii with increasing subhalo age. The mean masses M200 in- 
ferred from the GGLC signals also reflect the mass loss. 

The GGLC profiles also indicate that the mass loss is not due 
to a simple truncation but rather due to a mass loss at all radii (in 



^ Only for the subhalo samples with the largest infall mass range, the 
decrease in subhalo mass can be partially attributed to differences in the 
mean infall mass with time of the samples. 
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Age [Gys] 


Mi„f [10'«Me//i] 


Mm [Mo] 


age 


X — 


t<3 




4.621 X 10" 


1.797 


+ — 


3 < f < 4.5 


20 < Mm < 60 


4.711 X 10" 


3.723 


• — 


t > 4.5 




4.741 X 10" 


5.608 


X- - 


t<3 




1.278 X 10'2 


1.852 


+- - 


3 < f < 4.5 


60 < Mm < 150 


1.279 X 10'- 


3.720 


• - 


t > 4.5 




1.269 X 10'- 


5.629 


X • • • 


t<3 




1.28 X 10" 


1.798 


+ . . . 
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Mm > 150 


6.82 X 10'^ 


3.688 


«... 
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4.88 X 10'^ 
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Fig. 10: Profiles for different infall mass and time spend inside the clus- 
ter: with error bars the regions where the bias is below 10%. The color 
and the line type are the same for samples selected from the same in- 
fall mass range. The symbols are the same for samples defined from 
the same age bin. In solid blue we plot the average signal for the host 
cluster as a visual reference. In the table the ranges used for the binning 
and the average value for each sub-sample. 



the range where the subhalo profiles can be measured reliably), 
which corroborates the findings of'Hayashi et al. (2003). At very 
large radii, the subhalo profiles tend to increase with increasing 
subhalo age. This result indicates that the surrounding subhalo's 
contribution and main halo asymmetries become more important 
with increasing subhalo age. 

5.3. Forecasts for future surveys 
5.3.1. Galaxy observables 

Understanding how different observables depend on the under- 
lying subhalo characteristics is essential for obtaining optimal 
results with our method. We study these relations using the semi- 
analytic models of galaxy formation in the MS. A discussion of 
how well the model predicts such relations is beyond the scope 
of this work. Our aim is to show the feasibility of extracting in- 
formation with the method we propose. A good review about 
semi-analytic models is d iscussed in|Baugh|(|2006| l and [Benson] 
( 2010[ l. Gal axy clustering ([Blaizot et al. 2006)1, galaxy colors and 
metallicity ( Springel et al.j2001| and morphology evolution ( De 
|Lucia et al.|2006| l,"just to outHne a few examples, have been suc- 
cessfully studied with semi-analytic models. 



We only include type-1 galaxies (galaxies with a host sub- 



halo) in our study. We justify our choice in Sect. 6.2 

For our investigation we have studied proxies for the mass 
as determined by SUBFIND, the infall mass and the age of the 
subhalos (i.e. the time spent inside the cluster). The results are 
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Table 8: Parameters for the different surveys considered in our de- 
tectability study. We present the median redshift, the background galaxy 
number density, the standard deviation of the intrinsic ellipticity and its 
equivalent survey solid angle. 



shown in Fig. [TTj where we have marked with a black frame 
the proxies that show the strongest correlations and are there- 
fore used for our final analysis, for which we only used galaxies 
hosted by subhalos and imposed an apparent magnitude limit of 
r < 22. 

The age of the subhalos is the most difficult to obtain; none of 
the observable that we considered displays a strong coiTelation 
with this quantity. In top left panel of Fig. 1 1 we plot the pro- 
jected radial separation of the subhalo versus its age. Subhalos 
fall gradually into the main halo through dynamical friction, 
and so c/m-s should be an indicator of the age. However, since 
the projected separation corresponds to a wide range of three- 
dimensional separations, we only find a very weak correlation. 

In the top right and the middle left panels, we plot the r-i and 
u — r color indices as a function of subhalo age. The evolution 
of the stellar population and the depletion of gas to form new 
stars, make the color of galaxies evolve with time. We account 
for the redshift evolution of the galaxy colors as follows: we 
identify for each redshift the red cluster sequence averaged over 
all clusters, and subtract its mean color from all galaxy colors 
at the same redshift. By doing this, we can compare galaxies at 
different redshifts. Both observables are only weakly correlated 
with the subhalo age. Note that the contour levels are not equally 
spaced. 

The most sensitive proxy for the age of the subhalos is their 
morphology, which we quantify using the ratio between the lu- 
minosities of the bulge and the disk (middle right panel). There 
are many galaxies which only have a disk component, or only 
have a bulge (25% of the subhalos in total). In both cases the 
age of these subhalos cannot be determined. We exclude these 
galaxies from our analysis, as they should be easy to avoid in 
real surveys. 

The commonly used proxy for mass is luminosity. In the 
bottom left panel, we plot the r band absolute magnitude as a 
function of Msubf- Luminosity in the / bands give qualitatively 
similar results. 

Finally, using results from |Wang et al.| ( |2006[ ), we study the 
stellar mass as a proxy for the infall mass Minf. The result is 
displayed in the bottom light panel, showing a strong correla- 
tion. In the semi-analytic models, galaxies loose their gas rapidly 
once they become part of a cluster and the star formation ef- 
fectively stops. Therefore galaxies in clusters maintain on aver- 
age the stellar mass which they had prior to fall into the cluster, 
which is then strongly coiTelated with the halo mass at that time. 

5.3.2. Simulated Surveys 

We now investigate the prospects of measuring the weak lensing 
signal of subhalos in upcoming surveys. We focus in particular 
on DES and LSST, but we explore also variations of the main 
survey characteristics. The area simulated in our ray-tracing is 
not large enough to represent the total size of the surveys. We 
are not able, therefore, to reproduce the expected cosmic vari- 
ance. However, our major concern is the detectability of the sig- 
nal due to the statistical noise, since we expect that our subhalo 
samples to be representative. Therefore the influence in the mea- 
surements stemming from cosmic variance is assumed to be sub- 
dominant. 

We create mock galaxy catalogues similar to those expected 



from these surveys. As described in Sect. 4.2 the redshift distii- 
bution of our sources follows that of IBrainerd et al. Ill 99 6 The 
modulus of the intrinsic ellipticity is drawn from a Gaussian dis- 
tribution with zero mean and standard deviation a-\e^\. This in- 
trinsic ellipticity is the dominant source of noise in our mea- 



survey 


Median z 


gal 
arcmin^ 


kil 


eff. degrees^ 


LSST 


1.2 


40 


0.05 


18432 


INT 


0.9 


25 


0.1 


4608 


DES-WIDE 


0.68 


12 


0.05 


18432 


DES 


0.68 


12 


0.1 


4608 



surement, resulting in an eiTor in the measurement oc cr|^.|/ ^|n, 
where n is the total number of foreground-background galaxy 
pairs used to measure the final A2;(^). For a single galaxy, we 
assume a-\^^\ - 0.3. However, the survey areas of DES and LSST 
are considerably larger than the area covered by the ray-tracing. 
In order to simulate these surveys with their larger number of 
background galaxies, we therefore use a effective intrinsic ellip- 
ticity variance, which is reduced by the ratio of the survey area 

to the area of the simulations ') = cr^f.| x ^■ 

For the lens sample, we impose a maximum redshift of z = 
0.9. Subhalos at higher redshift are difficult to probe with lensing 
since there are fewer background sources available. Henceforth, 
their inclusion offers no improvement in the detectability of the 
signal. For the LSST-like survey, we only considered lens galax- 
ies with r < 26 in apparent magnitude, for the other three surveys 
is r < 22. We divide these samples further into 14 bins in abso- 
lute r-band magnitude ranging from r - -23.66 to r = -19 with 
a bin width of 0.33 magnitudes. 

The details of these surveys are summarized in Tab. [8] The 
LSST and the DES surveys are expected to be around 10% larger 
in size than the ones that we effectively consider here. In this 
way, we account for the reduction of the usable area due to mask- 
ing. Furthermore, with the ficticious DES-WIDE survey we an- 
alyze the impact of having a survey like DES but with the solid 
angle of LSST, and with the INT survey the effect of having 
deeper observations on the same area as DES. 

For simplicity, in the following the errors plotted are the 
combined standard errors obtained from Ssubh and Scai. We as- 
sessed that the correlation between Esubh and Scat is subdomi- 
nant. However, for the parameter fitting we use the full covari- 
ance matrix computed with bootstrapping in each case. 



5.3.3. Signal-to-noise ratio 

For each survey and luminosity bin, we compute a mean signal- 
to-noise ratio (SNR) by averaging over eight logarithmic bins in 
projected radius ^ in the range from 0.043 to 0.12 Mpc: 



SNR = 



gt(^) 



(32) 



We use for cr [et(^)] the combined standard eiTor of the mean 
2^Subh and Seal. These SNR values neither take into account cor- 
related noise nor express how well a particular feature can be 
detected. Nonetheless, they are compact estimators and offer a 
rule of thumb for the expected signal quality. 



In Fig. 12 we present the resulting SNR values for the four 
surveys considered. Luminous galaxies usually have a larger 
halo with a stronger signal, whereas faint galaxies are more nu- 
merous, corresponding to a smaller statistical error. The SNR 
depends on these two competing effects, resulting in a maximum 
of the SNR signal-to-noise ratios are around -21.5 magnitude. 
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The least luminous bins are below the two-sigma level in all 
surveys except for LSST; in these cases, the signal-to-noise ratio 
can be improved by using broader bins. 



Finally, in the left panel of Fig. 13 we present the predicted 
galaxy-galaxy lensing signals for a LSST-like survey, and in the 
right panel for a DES-like survey. The signal-to-noise is lower 
in a DES-like survey, requiring the use of fewer luminosity bins 
(from r - -24.25 to r = 19 with a width of 0.75 magnitudes) 
than for LSST. This figure presents a global characterization of 
the expected signals, taking into account an adequate treatment 
of the galaxy samples. We can see the quality of the signals and 
how detailed our models can be. 



5.3.4. Mass-luminosity relation 

We now address the question how well a mass-luminosity re- 
lation for subhalos can be estimated from DBS and LSST. The 
measurement of such a relation can help to distinguish between 
models of sub-structure formation. As we probe larger scales 
than stellar dynamics, and we consider large samples of galax- 
ies compared to strong lensing analyses, our method is largely 
complementary to other methods. In Fig.[T3]we present the re- 
sults for DBS and LSST. Note that we only plot error bars for 
the range where we consider our measure has a bias below 10%. 

In order to obtain the subhalo mass Msubf, one needs to de- 
rive M200 from the measured AE profiles. The subhalo mass can 
then be computed from a scaling relation between Msubf and 
M200, which is shown in the top panel of Fig. 14 The final mass- 
luminosity relations, for both a LSST and a DBS-like survey, are 
shown in Fig. 15 Using DBS or LSST, one can constrain the 



mass-luminosity relation over two decades in mass. The system- 
atic differences between the relation for the two surveys can be 
explained by the different lens samples used. The brightest bins 
are populated with the same objects for both surveys and there- 
fore the scaling relations are similar At the fainter end, the dif- 
ferent magnitude cuts in combination with the broader luminos- 
ity bins in the DBS change the mass-luminosity relation (Fig. 14 
bottom panel). 

Nonetheless, the relation between M2()o and the mass is the 
same (within the errors) in both cases. The smaller number den- 
sity of faint galaxies in DBS compared to the LSST survey re- 
duces in this case the range over which we can measure the mass- 
luminosity relation. Note that the faintest luminosity bin for the 
DBS-like survey may not be usable as the measurement error 
becomes too large. 

5.3.5. Evolution of subhalo profiles 



In Sect. 5.2.3| we showed how the subhalo mass profiles change 
with time and found that it is not possible to detect a truncation 
in the original sense, and that the mass density decreased at all 
scales consistent with a heating process. We now study whether 
this behavior can be detected with DBS or LSST using observ- 
able quantities only. 



In Fig. 16 the predicted galaxy-galaxy lensing signals for 
both LSST (left panel) and DBS (right panel) are shown with dif- 
ferent symbols for the different stellar mass bins. Additionally, 
we have divided each of the stellar mass bins into two subsets 
according to morphology. The details of the binning in morphol- 
ogy and stellar mass are given in tables inside the figure. 

The dashed red lines belong to galaxies with a relatively 
large bulge, which should belong to old subhalos. With blue dot- 
ted lines we plotted the signals from galaxies with a large disk. 



which should belong to younger subhalos. The number of stel- 
lar mass bins we can analyze is limited, since we need to split 
them further into early-type and late-type galaxies. The stellar 
mass ranges were obtained by experimentation; in the case of 
the DBS-like survey we focused on a particular range where the 
effect was strongest. Note that we only plot error bars for the 
range where our measurement has a bias below 10%. 

The sub-division according to morphology does not always 
produce the desired results. One can see from the last two rows 
in the top table, that for the most massive satellite galaxies, the 
infall mass (eighth column) of the two morphology bins are very 
different. Older subhalos have a significantly larger mean infall 
mass, which indicates that they evolved in a different manner. 
This difference masks any possible detection of mass loss due to 
tidal stripping by the cluster. 



With the exception of this particular bin, however. Fig. 16 
shows that the profile amplitude of bulge-dominated galaxies is 
smaller than the one of disk-dominated objects, in accordance 
with our predictions from Sect. 5.2.3 This can be detected with 



high significance in the LSST; for the DBS-like survey the de- 
tection is hampered by the low signal-to-noise ratio. 



6. Systematic effects 

6.1. Validity of tlie weak lensing approximation 

So far, we have assumed that the weak lensing approximation is 
valid in our measurements, i.e. that <k 1 and thus that the re- 
duced shear effectively equals the shear. To check the validity of 
this approximation, we compare the galaxy-galaxy lensing sig- 
nal of subhalos measured using shear with the signal obtained 
using reduced shear. 

In Fig. [14] we have estimated M200 by fitting NFW profiles 
to the AL profiles measured from galaxy catalogues obtained 
from the full ray-tracing using the shear yiO) and from similar 
catalogues using the reduced shear g{0). No significant differ- 
ence between the two methods is found. For the average subhalo, 
the impact of the environment is not important. Within internal 
checks, we have observed a significant discrepancy only for the 
most massive subhalos for the innermost radii. This is produced 
by the large convergence of the subhalo itself and not related to 
the host main halo. 



6.2. Type-2 galaxies 



The semi-analytic model of galaxy formation by De Lucia & 
[Blaizot (2007 ) contains galaxies that have "lost" their subhalos. 
Whenever a subhalo is no longer detected (e.g. because it has lost 
most of its mass due to tidal stripping), the hosted galaxy is still 
maintained and placed at the position of the most bound particle 
of the former subhalo. After the dynamical friction timescale, 
these galaxies eventually disappear as they merge into the main 
halo galaxy. 

The results presented so far were computed excluding type-2 
galaxies for the following reasons: 

- they are hosted by subhalos below or around the mass reso- 
lution limit, 

- they exist under the assumption that the stripping of a large 
part of the subhalo does not destroy the hosted galaxy, 

- the survival time is heavily influenced by the model. 

Nevertheless, to understand the impact of such galaxies on 
our measurements, we analyze how our results are altered when 
including them. We present how our measurements would be 
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affected assuming that type-2 are real and that we cannot distin- 
guish between them and type-1 galaxies. 

In Fig. 17 we present AS(^) for a selection of luminosity 
bins for type-1 galaxies (with dashed blue lines) and the same 
luminosity bins but for the combination of type-2 and type-1 
galaxies (with solid black lines). The inclusion of type-2 galaxies 
change the average profile, especially at large radii. The change 
is larger for fainter galaxies. The new profiles are flatter, and 
there is a systematic error in the determination of the mass pro- 
file. This effect can also interfere with our ability to assign any 
kind of mass estimate. 



In Fig. 18 we present the same analysis as in the top panel 
in Fig. 16 but with both type-2 and type-1 galaxies. Recall that 
the line and color type encodes the morphology (as a proxy for 
the subhalos' age), and the symbol encodes the stellar mass (as a 
proxy for the subhalo infall mass). The inclusion of type-2 galax- 
ies is more important if we aim to estimate the mass loss on sub- 
halos induced by tidal stripping. Our former ability to measure 
this effect was small, and the change in the profiles that type-2 
galaxies produce makes it impossible. 

We are able to trace the mass responsible for the lensing sig- 
nal in type-2 galaxies. In Fig. 19 we quantify how likely is to 



find a companion near a type-2 galaxy. We compare the number 
density as function of radius, around a galaxy and its calibration 
point. Type-2 galaxies are on average close to a subhalo which 
accounts for the lensing signal measured. We were able as well 
to detect type-2 galaxies associated with projected overdensities 
which were not detected as subhalos by SUBFIND. 

The presence of galaxies with no halos makes our measure- 
ments less precise. However, if such galaxies exists, we could 
use our analysis to determine the fraction of galaxies of such 
type, and to distinguish them from galaxies with a host halo. 

7. Summary and conclusions 

In order to understand the evolution of galaxy clusters, it 
is essential to study the matter profiles of satellite galaxies. 
Analyzing how these profiles evolve with time can improve our 
knowledge on galaxy evolution and on structure formation. 

In this work we analyzed the use of weak gravitational lens- 
ing on satellite galaxies inside clusters, in particular the use of 
galaxy-galaxy lensing. With this probe, we can measure in a non- 
parametric way the average projected mass profiles of the host 
matter halos of galaxies. This cosmological probe correlates the 
image distortion (shear) of a background galaxy with the mass 
of a foreground galaxy. 

Galaxy-galaxy lensing needs large galaxy samples. For this 
reason it has not yet been fully exploited on cluster satellite 
galaxies. In our work, we forecast results for future surveys us- 
ing the Millennium Simulation ( Springel et al.|2005] l, ray-tracing 
simulations (Hilb ert et al.||2009 l and the galaxy semi-analytic 
catalogues by De L ucia & Blaizot| (| 2007| . 

We presented the details of galaxy-galaxy lensing measure- 
ments on satellite galaxies and its contamination from the host 
cluster. In order to overcome the contamination from the host 
cluster we proposed a calibration method, which can solve the 
problem up to a certain range in radius. For each subhalo, we 
define a point at the same distance from the main halo as the sub- 
halo, but in the opposite direction as seen from the cluster center 
We can estimate the contribution of the main halo by measuring 
the tangential shear around this new point, under the assumption 
that the cluster is point-symmetric. 

We created mass maps of our clusters using theoretical pro- 
files with known parameters. With this mock cluster sample we 



were able to test our calibration method taking into account re- 
alistic characteristics for our cluster samples such as halo spatial 
distribution or mass function. We estimated the performance of 
our measurements at different radii and we defined a minimal 
separation between the subhalo and the main halo center to op- 
timize the signals. 

With the previous tests, we could characterize the subha- 
los in the Millennium Simulation using projected mass maps. 
The weak lensing signal of subhalos is well described by a sim- 
ple NFW profile, and it was not possible to detect a truncation 
radius. Our results are consistent with an abrupt truncation of 
the mass profile at radii larger than 0.2 Mpc. There are certain 
discrepancies between our work and the previously published 

(|2007l). 



works from Limousin et al. 



Natarajan et al.| ( |2007| l. 



Halkola et al.| ( |2007| l and |Suyu & Haikola| ( |2010| l. These authors 



measured a much smaller extent of the subhalos using gravita- 
tional lensing on a few observed clusters. The results were de- 
rived using parametric models for the mass profiles of the subha- 
los and the main halo. The models used by these authors do not 
fit subhalos in the Millennium Simulation. Since the Millennium 
Simulation only contains dark matter, there is not absolute cer- 
tainty that NFW profiles should be the best description for ha- 
los of real galaxies. On the other hand our method did not in- 
clude strong lensing constraints in comparison to the previous 
works. Nevertheless, our results challenge the choices made in 
the afore-mentioned works. The truncation radii that they mea- 
sured can be also interpreted as a consequence of the paramet- 
ric method used. Therefore, further analysis and better data are 
needed in order to solve the problem. 

We also were able to characterize the evolution of the mass 
profiles. We show that the lensing profiles decrease in amplitude 
with time. This is consistent with a mass loss at all scales. Due 
to the tidal forces exerted by the host halo, subhalos are stripped 
of the mass at the outermost radii, and at the same time the mass 
at the inner regions is redistributed. This is consistent with the 



work by Hayashi et al. ( 2003, ) and also supports the idea of large 
truncation radii. 

After describing the subhalo profiles we used simulated 
galaxy catalogues to forecast signals for future surveys, focus- 
ing on PES (|The D ark Ene rgy Survey Collaboration|2005[ ) and 
LSST ( |Ivezic et~al. 1,200 8 ). We analyzed the semi-analytic cat- 
alogues in order to classify the galaxy samples to optimize the 
measurements. With the result from this analysis we derived the 
following results. 

We predict the detectability of the signals using a compact 
estimator of the signal-to-noise ratio. The cluster sample re- 
quired is very large, but we can already expect signals from a 
DES-like survey roughly above the three sigma level. The data 
from a LSST-like survey is well-suited for the studies that we 
proposed. 

There is not a unique way of separating between subhalo and 
host halo mass. In this work, we considered that the mass of the 
subhalo is well estimated in our simulations by the SUBFIND 
algorithm ( [Springel et al.|200l) (Msubf)- We modeled the rela- 
tion between the measured NFW profiles and the gravitationally 
bound mass of the subhalo mass Msubf- We also checked that 
the weak lensing approximation is valid. 

According to the semi-analytic catalogues, luminosity in the 
SDSS r band is a good proxy for mass. We predict that it is possi- 
ble even for a DES-like survey to constrain the mass-luminosity 
relations of subhalos over two decades in mass, from around 
5 X 10" Mq to lO'-* Mq. 

In order to study the time evolution of the profiles using re- 
alistic data, we binned the galaxies according the mass of the 
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subhalo prior to falling into the cluster and the time spent inside 
the cluster (subhalo age). Within the semi-analytic catalogues, it 
is possible to infer their initial mass from their observed stellar 
mass. On the other hand we could not find any galaxy observ- 
able that is strongly coiTelated with the subhalo age. With these 
observables, we could only put weak constraints to the evolution 
of subhalo matter profiles. 

In our analyses we neglected galaxies without a host subhalo 
which one finds in the semi-analytic catalogues of De Lucia & 
Blaizot (2007 1. A priori, these galaxies are not fully reliable as 
they populate mass overdensities below the resolution limit of 
the Millennium Simulation. Neveitheless, for completeness we 
considered their influence. These galaxies show a lensing signal 
with a high amplitude, which can be explained as being produced 
by correlated halos. We also were able to quantify the correla- 
tion between the position of type-2 galaxies and other subhalos. 
Finally, we presented how our previous analysis is affected by 
assuming that these galaxies are realistic and that we are not able 
to distinguish them from galaxies with a host subhalo. However, 
further simulations are needed in order to investigate galaxies of 
such characteristics. 
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Fig. 11: Relations between physical subhalo properties (vertical axis) and observable proxies (horizontal axis). With black frames, we outline 
those finally used. Contours give the full distribution; we overplot the mean value of each subhalo property and its standard deviation for each bin 
in the observable (curves with error bars). Top left: age of the subhalo versus the projected radial separation of the subhalo from the main halo 
center. Top right: age of the subhalo versus redshift-corrected r - i. Middle left: age of the subhalo versus redshift-corrected ii - r. See text for a 
description of the redshift correction. Middle right: age of the subhalo versus our morphology estimator. Bottom left: Msubf versus absolute r band 
magnitude. Bottom right: versus stellar mass (using results fromlWang et al. 2006 1. Note that the contour levels differ between the panels. 
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Fig. 12: Our compact SNR estimator, for the different luminosity bins 
for each mock survey considered. For each luminosity bin we plot the 
mean magnitude. The errors plotted are standard errors obtained from 
the sample. 
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Physical radius [Mpc] Pliysical radius [Mpc] 

Fig. 13: Measured AX(^) profiles of galaxies binned by r-band absolute magnitude. Left panel: LSST-like survey, with r < 26 in apparent 
magnitude. Right panel: DES-like survey, with r < 22 in apparent magnitude. Error bars show the data used for the analysis later on, the range 
was derive from our results in Tab. [51 The solid blue line without error bars is the measurement for the host cluster shown as a visual reference. 
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Fig. 15: Bottom: the relation between the fitted parameter M200 and 
the luminosity in absolute observer frame r band magnitude. We plot 
in solid black the results for a DES-like survey, and in dashed red for a 
LSST-like survey. 
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Fig. 14: Top: relation between the fitted parameter M200 and Msubf- 
Bottom: the mass-luminosity relation obtained from the catalogues. We 
plot in solid black the results for a DES-like survey, and in dashed red 
for a LSST-like survey. 
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Data for the LSST-like survey (left panel). Galaxies with z < 0.9 and r < 26. Masses in units of 10'" Mg 
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Data for the DES-like survey (right panel). Galaxies with z < 0.9 and r < 22. Masses in units of 10'" M( 
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Fig. 16: Morphology-stellar mass classification for a LSST-like survey (left panel) and for a DES-like survey (right panel). The black solid line is 
the measurement for the host cluster shown for visual reference. We only plot symbols and error bars for the range where we consider that the bias is 
below 10%. The subhalos are at ^m-s < 0.5 Mpc. Blue dotted lines correspond to galaxies with a relatively small bulge: Disk — > < Lsuige/i'Totai ^ 
0.6. Red dashed lines correspond to galaxies with a large bulge which spent more time inside a cluster: Bulge — > 0.6 < LBuigc/i'Totai < 0.98. The 
different symbols distinguish different stellar mass bins (in units of 10'" Mq). For the left panel the ranges are: 1 — » [0.14 : 1.4], II — > [1.4 : 4.1 1], 
IIl-» [4.11 : 10.96], 4^ [10.96 : 00]. For the right panel the ranges are I -> [1.1 : 4.1], 11^ [4.1 : 13.7]. 
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Fig. 17: Excess surface mass density AE(^) for a few luminosity bins 
for an LSST-like survey. The blue dashed lines correspond to the mea- 
surements only for type-1 galaxies. The solid black lines correspond to 
the combination of both type-1 and type-2 galaxies. The range where 
we assume the measurement has a bias below 10% is highlighted, (de- 
rived using our results for type-1 only profiles). The solid thick red line 
is the signal for the average cluster, shown for visual reference. 
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Fig. 18: Morphology-stellar mass classification for a LSST-like survey 
including both type-1 and type-2 galaxies. Blue dotted lines correspond 
to galaxies with a relatively small bulge, red dashed to those with a large 
bulge. The different symbols distinguish different stellar mass bins. The 
range where we assume the measurement has a bias below 10% is high- 
lighted, (derived using our results from the type-1 only profiles). The 
classification ranges are described in the top table in Fig.[T6] The black 
solid line is the measurement for the host cluster shown for visual ref- 
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Fig. 19: Ratio between the galaxy number density around a galaxy and 
around its corresponding calibration point, as a function of distance. We 
also plot as visual reference the line where the ratio is 1 . Dotted-dashed 
red: excess type-2 number density around type-2. Solid blue: type-1 
number density around type-1. Green dashed: excess type-1 number 
density around type-2. ^m-s < 0.5 Mpc. 



